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Аннотация 

Введение. Для всех полимерных материалов и композитов на их основе характерны явно выраженные 
реологические свойства, прогнозирование которых является одной из важнейших задач механики полимеров. 
Большие возможности для прогнозирования реологических параметров полимеров открывают методы 
машинного обучения. Ранее проводились исследования на предмет построения прогнозных моделей с 
использованием искусственных нейронных сетей и алгоритма СаВооз(. Наряду с этими методами, благодаря 
возможности обрабатывать данные с сильно нелинейными зависимостями между признаками, широкое 
применение в смежных областях находят методы машинного обучения — метод А-ближайших соседей и метод 
опорных векторов (5УМ). Однако ранее к проблеме, рассмотренной в данной статье, эти методы не применялись. 
Целью работы явилась разработка прогнозной модели для оценки реологических параметров полимеров 
методами искусственного интеллекта на примере поливинилхлорида. 

Материалы и методы. В работе применены метод А-ближайших соседей и метод опорных векторов для 
определения реологических параметров полимеров на основе кривых релаксации напряжений. Обучение 
моделей выполнялось на синтетических данных, сгенерированных на основе теоретических кривых релаксации, 
построенных с использованием нелинейного уравнения Максвелла-Гуревича. Входными параметрами моделей 
выступали величина деформации, при которой производился эксперимент, начальное напряжение, напряжение 
в конце процесса релаксации, время релаксации и условное время окончания процесса. Выходные параметры: 
модуль скорости и коэффициент начальной релаксационной вязкости. Модели разработаны в среде 
7аруег МоеБоок на языке Ру®оп. 

Результаты исследования. Построены новые прогнозные модели для определения реологических параметров 
полимеров на основе методов искусственного интеллекта. Предложенные модели обеспечивают высокое 
качество прогнозирования. Метрики качества модели в алгоритме 5УВ составляют: МАЕ — 1,67 и 0,72; 
МЕ — 5,75 и 1,21; ВМЗЕ — 1,67 и 1,1; МАРЕ — 8,92 и 7,3 для параметров начальной релаксационной вязкости 
и модуля скорости соответственно с коэффициентом детерминации А? — 0,98. Разработанные модели показали 
среднюю абсолютную процентную ошибку в диапазоне 5,9-8,9 %. Помимо синтетических данных, 
разработанные модели также апробировалась на реальных экспериментальных данных для поливинилхлорида в 
диапазоне температур от 20° до 60 °С. 

Обсуждение и заключение. Апробация разработанных моделей на реальных экспериментальных кривых 
показала высокое качество их аппроксимации, сопоставимое с другими методами. Таким образом, алгоритмы 
К-ближайших соседей и ЗУМ могут использоваться для прогнозирования реологических параметров полимеров 
как альтернатива искусственным нейронным сетям и алгоритму Са Вооз$ь требующая меньших усилий по 
предварительной настройке. При этом в данном исследовании наиболее предпочтительным методом машинного 
обучения оказался метод ЗУМ, так как он более эффективен в обработке большого числа признаков. 
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опорная векторная регрессия 
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Введение. Полимеры применяются в различных областях промышленности, включая производство пластика, 
текстиля, упаковочных материалов и многих других. Точное прогнозирование реологических параметров 
полимеров является сложной задачей, имеющей важное значение для оптимизации производственных процессов 
и создания продуктов с желаемыми свойствами. 

На сегодняшний день методы машинного обучения приобрели большую популярность в различных областях, 
включая химию и материаловедение, благодаря своей способности эффективно обрабатывать и анализировать 
большие объемы данных. Данные методы позволяют прогнозировать свойства материалов. В работе [1] описана 
платформа, основанная на машинном обучении, предложена интеграция метрологического обеспечения в 


Механика 


Б&р://уезштК-Чопзва.га 


'А4уапсеа Епятеетия Везеагсй (Козоу-оп-Оопн). 2024;24(1):36-—47. е155М 2687—1653 


условиях цифровой трансформации. В [2] прогнозируется локальное распределение деформации, развитие 
пластической анизотропии и разрушения в аддитивно изготовленных сплавах. Проблемы разработки регуляторов 
измерительного контроля на цифровых платформах сформулированы в работе [3]. Интеллектуальная модель 
управления параметрами сварки соединений внахлестку построена в [4]. Однако вопросы применения методов 
машинного обучения для прогнозирования реологических свойств полимеров остаются недостаточно 
исследованными. Это вызвано как техническими, так и методологическими сложностями, такими как 
неоднородность структуры полимеров, их чувствительность к внешним условиям и сложные взаимодействия 
между молекулами в процессе деформации. 

Исследования в области реологических свойств полимеров и композитов методами машинного обучения 
несут большие перспективы в строительной индустрии [5]. Для многих полимеров экспериментальные данные 
хорошо описывает обобщенное нелинейное уравнение Максвелла-Гуревича [6], для одноосного напряженного 
состояния имеющее вид [7]: 


д5* _ Л 
0 п“ 
Г* =в- ЕЕ", (1) 
Е а 
-  . @Р| „| 
т 1 т 
где 5’ — деформация ползучести; /’ — функция напряжений; с — напряжение; Е» — модуль 


высокоэластичности; Ио — начальная релаксационная вязкость; т” — модуль скорости. 

Для определения реологических параметров полимеров, таких как начальная релаксационная вязкость (далее 
просто «вязкость») и модуль скорости, можно использовать различные интеллектуальные модели машинного 
обучения [8, 9]. Например, одной из таких моделей является нейронная сеть, которая может обучаться на 
сгенерированных наборах данных для определения оптимальных параметров полимеров [10]. 

Прогнозирование на синтезированных данных достаточно распространённая практика, в том числе для 
методов нелинейной оптимизации [11, 12]. Одним из способов генерации данных выступает применение 
функций Розенброка, Химмельблау, Бута [13], которые эффективно используются для тестирования методов 
оптимизации, таких как методы градиентного спуска, генетические алгоритмы и метод Ньютона. Такой подход 
применен в работе [14], где для проверки эффективности различных методов оптимизации сгенерирован набор 
данных на основе теоретических кривых релаксации напряжений с использованием нелинейного уравнения 
Максвелла-Гуревича. 

В работе [15] для прогнозирования долговечности железобетонной балки приведено несколько подходов 
машинного обучения, такие как нейронная сеть обратного распространения, линейная и гребневая регрессия, 
дерево решений, случайный лес. Входными параметрами исследования явились как различные характеристики 
материала, так и их свойства, в зависимости от окружающей среды (температура, влажность). В итоге по 
результатам исследования модель обратного распространения определила более точный прогноз (85 %), средние 
значения (МАЕ) и МАРЕ составили 1,13 % и 14,5 % соответственно. 

Другой подход к решению обратных задач теории ползучести методом нейронных сетей основан на обучении 
модели на больших объемах экспериментальных данных. В работе [16] была разработана нейросетевая модель, 
которая обучалась на данных, полученных в результате многолетних экспериментов над полимерными 
материалами, и успешно прогнозировала вязкоупругое поведение этих материалов. Для исследования были 
использованы данные, полученные в результате экспериментов над образцами различных материалов. 

В отличие от вышеперечисленных работ, представленная работа призвана содействовать развитию более 
точных и надежных методов прогнозирования свойств полимеров, таких как метод А-ближайших соседей и метод 
опорных векторов, что имеет важное значение для различных отраслей промышленности и науки. 

Цель исследования заключалась в разработке прогностической модели на основе методов искусственного 
интеллекта для анализа реологических свойств полимеров. Ранее авторы в своих работах уже применяли 
алгоритм машинного обучения на основе градиентного бустинга СаВооз( для обработки кривых релаксации 
напряжений [17, 18]. Са Воозё является одним из наиболее мощных алгоритмов машинного обучения, 
применимым для решения не только задач регрессии, но и задач классификации и ранжирования [19]. 

Метод Са{Вооз{ может быть полезным для решения некоторых задач, но он также имеет свои ограничения и 
недостатки. В связи с этим появился интерес к использованию и других, упомянутых ранее [20], алгоритмов для 
решения поставленной задачи. 
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Материалы и методы. Сгенерированный массив данных частично представлен в таблице 1. Данный массив 
был сформирован на основе теоретических кривых релаксации напряжений, описываемых уравнением 
Максвелла-Гуревича, по методике, представленной в работе [14]. Диапазоны изменения модуля скорости и 
начальной релаксационной вязкости в сгенерированном массиве соответствуют реальным диапазонам для 
поливинилхлорида в интервале температур от 20° до 60 °С. Общее количество численных экспериментов (и) 
составляло 30 000. 


Таблица 1 
Таблица исходных данных для обучения модели 
Напряжение | Напряжение ВЕ Условное ме Вязкость 
№ ее дих, в начале в конце о акСадий время короетИ п&, 10° 
% процесса со, | процесса о», р окончания т”, МПа а 
МПа МПа процесса 195, ч ы 
1 1,000 190,000 0,909 0,277 1,484 6,000 3,000 
2 2,000 20,000 1,818 0,109 1,003 6,000 3,000 
3 3,000 30,000 2,727 0,046 0,820 6,000 3,000 
4 1,000 10,000 0,909 0,861 4,615 6,000 9,333 
5. 2,000 20,000 1,818 0,339 3,122 6,000 9,333 
6 3,000 30,000 2,727 0,142 2,552, 6,000 9,333 
7 1,000 190,000 0,909 1,445 7,747 6,000 15,667 
29997 3 45 37,5 0,285 2,476 15 53,666 
29998 1 15 12,5 1,003 4,255 15 60 
29999 2 30 25,0 0,558 3,371 15 60 
30000 3 45 37,5 0,319 2,769 15 60 


Набор данных составил пять входных переменных и две выходных переменных. Входные переменные 
(единица измерения): деформация — = (%); напряжение в начале процесса — оо (МПа); напряжение в конце 
процесса — о» (МПа); время релаксации — (1, (ч); условное время окончания процесса — 155 (ч). Выходные 
переменные (единица измерения): модуль скорости — т” (МПа); начальная релаксационная вязкость — По 
(в таблице 1 и далее просто «вязкость») (10° МПа-с). Величины о0, ох {, и 65 схематически показаны на типовой 
кривой релаксации напряжений, представленной на рис. 1. 


Рис. 1. Типовая кривая релаксации напряжений 


Алгоритм А-ближайших соседей (А-ММ) основан на анализе сходства ближайших объектов. Метод А-ММ 
пользуется большим спросом при решении различных типов задач машинного обучения. 
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Формула (2) представляет общий вид алгоритма, где и(р, х) — весовая функция, оценивающая степень 
важности 1-го соседа. 


Е(х)=агетахусу Хы =У] и(ьх), _ 


Максимально суммарный вес может достигаться для нескольких объектов одновременно. Энтропию этого 
процесса можно регулировать с помощью нелинейной последовательности и(ьх)=[1< а’ (метод 
экспоненциально взвешенных А-ближайших соседей) при условии, что 0 < 4 < 0,5. 

Представляя собой достаточно простой алгоритм машинного обучения, К-ММ хорошо применим при решении 
задач классификации и регрессии. Преимущества данного метода: простота реализации, отсутствие необходимости 
в предварительном обучении модели. Он используется для любых типов данных, включая категориальные и 
числовые. Недостатки: склонность к переобучению (при условии, если К слишком мало), низкая 
производительность при больших объемах данных, нет возможности учитывать взаимосвязь между признаками. 

Алгоритм опорных векторов (ЗУВ) — регрессия опорных векторов, решает задачи минимизации суммы 
средней абсолютной ошибки. ЗУК является более устойчивым к выбросам, в отличии от метода наименьших 
квадратов, за счёт коэффициента регуляризации (С) и «эпсилон-нечувствительной трубки» (=). При этом & 
определяет ширину трубки, в которой ошибки игнорируются. Для нахождения минимума функции используется 
стохастический градиентный спуск. 

Алгоритм обучения по методу опорных векторов представляет собой функцию Р(х) аппроксимации и 
регуляризации эмпирического риска, преобразующую обучающую и тестовую выборки в выходные данные для 
каждого объекта соответствующей выборки. Формула (3) представляет общий вид алгоритма; (4) — линейно 
разделимая выборка; (5) — линейно неразделимая выборка, где С — коэффициент регуляризации; 
МХ», \о) — скалярное произведение векторов (признака и опорного вектора); и’; — весовые коэффициенты. 

Е(х)-сУ" (1-м, (ож) + -> тт Е(х). (3) 


У, Ио 

1 |2 . ` 

М > тие (х); 4) 
М (®т)>1={1:Д. 


МЕ +сУ и эттЕ(х); 


У, мо 
М; (*,т)>1-2,1={1:0; (5) 
2; >0,1={1:0. 


Функция К (8%) есть функция от пары объектов [9 представима в виде скалярного произведения в 
некотором пространстве Н, для которого имеет место преобразование \:Х-—Н. Функция 
К:ХхХх > В — ядро, если К (х,х') = (ч(х) м (х')), при условии, что К симметричная: К (х,х') = К (х,х') и 


неотрицательно определённая: [к [6х ) |4 (х) |4 (х') Чхах', \е:Х > К. Коэффициент  регуляризации 
определяется методом скользящего контроля. 

Преимущества метода ЗУМ: высокая точность в задачах классификации в нелинейных пространствах; 
способность работать с болышим количеством признаков (включая категориальные и числовые), обобщать 
данные (что позволяет применять модель для новых данных), работать с данными, которые не являются линейно 
разделимыми благодаря использованию ядерных функций. 

Недостатки метода ЗУМ: неэффективность работы с большими объемами данных; низкая 
интерпретируемость модели; требование настройки множества параметров, таких как тип ядра (его параметры, 
параметры регуляризации) и т.д. 

В данном исследовании алгоритмы разрабатываются в интеллектуальной вычислительной среде 
7аруег МоеБооК методами машинного обучения. 


т 


В качестве алгоритма обучения рассматривалась функция Ё(х), преобразующая обучающую {х,} ,ЕХ "и 
тестовую выборки т } ый еХ' в выходные данные при обучении { У ЕХ” и при тестировании { У у Ех" 


1= 11 
для каждого объекта соответствующей выборки. Обучение вектора параметров и; Е Й’ встроено внутри алгоритма. 


= 


В условиях представленной задачи: {у} = фактические значения вязкости 70 (в начале релаксационного 


т 
11 


т 
=1 


* 
процесса) и модуля скорости т”; { У }, — прогнозируемые значения вязкости По (в начале релаксационного 


процесса) и модуля скорости т”. 
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Подбор такого параметра, как количество соседей, влияет на обобщающую способность разработанной 
модели и является важным для еб корректной работы. Наиболее подходящий алгоритм расчёта расстояния на 
основе данных является ПО1$апсе, при котором веса объектов обратно пропорциональны их расстоянию. 
Соответственно, в случае более близких соседей объекта запроса, они имеют большее влияние, нежели их соседи, 
находящиеся на большем расстоянии от объекта. 

Набор данных разбит на обучающую и тестовую выборки в соотношении 75/25. В свою очередь, 20% 
обучающей выборки становится валидационной. Размерность выборок составила: обучающая — хина = 20 400; 
тестовая — хь»х = 6 000; валидационная — хе! = 3 600. Для переменных Урат, Уезь Уеми Данные распределены 
аналогичным образом. 

Для построения модели А-ближайших соседей были выбраны следующие параметры: количество соседей, 
размер листа, интервал, функция веса. Диапазон и функционал значений для настраиваемых параметров 
представлен в таблице 2. 


Таблица 2 
Таблица параметров для модели А-ММ 
№ Параметр Значение Функционал 
о Определяет оптимальное количество соседей 
1 Количество соседей (А) 3,5, 7,9 р 
для запроса 
Определяет скорость построения запроса 
2 Размер листа (и) 15, 20, 30 а о = . 
и необходимую память для хранения дерева 
Определяет параметр мощности (метрика 
3 Интервал (р) 1 (11), 2 (2) ге о т 
Минковского) 
4 Функция веса (у) ‘поп’, 'Ч5вапсе"' Прогнозирование весов 


Для построения модели ЗУК были выбраны следующие параметры: тип ядра, порядок ядра, коэффициент 
регуляризации (квадратичный регуляризатор), =. Диапазон и функционал значений для настраиваемых 
параметров представлены в таблице 3. 


Таблица 3 
Таблица параметров для модели ЗУК. 
№ Параметр Значение Функционал 
"Ппеаг'; 'ройу'; ТОР; Определяет тип гиперплоскости 
1 Тип ядра о и ев. й р 
9121014 (линейная/нелинейная) 

2 Порядок ядра 1,2, 3, 4, 5,7 Определяет степень полиномиальной функции ядра 
3 Квадратичный 2.3: 4:5. 7: 10 Решает проблемы мультиколлинеальности вектор 

регуляризатор (С) о и переобучения модели 
4 8 0,1; 0,2; 0,5; 1; 1,5;2; 3 Определяет отклонение объекта (мера близости) 


Результаты исследования 

На рис. 2 представлены корреляционные связи между переменными. 

Можно отметить следующие виды линейных корреляционных связей между отдельными входными и 
выходными переменными модели: 

— достаточно сильные — между переменными «Деформация» и «Напряжение в начале» рог = 0,93; «Время 
релаксации» и «Время окончания процесса» ри» = 0,93; 

— средние — между переменными «Деформация» и «Напряжение в конце» ро.. = 0,71; «Напряжение в 
начале» и «Напряжение в конце» рос. = 0,75; 

— слабые — между переменными «Время окончания процесса» и «Вязкость» р..,, = 0,58; «Вязкость» и 
«Время релаксации» р„., = 0,46. 

Наличие умеренной корреляционной связи между переменными или её отсутствие говорит лишь об 
отсутствии линейной связи, следовательно, возможно наличие нелинейной связи между переменными. 
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Рис. 2. Корреляционная матрица 


В таблице 4 предоставлены статистические характеристики исходного набора данных. 


Таблица 4 
Статистические характеристики исходного набора данных 
Параметр 5 60 С» И 195 т’ "о 

Ед. измерения % МПа МПа Ч Ч МПа 106 МПас 
соит 30 000,00 | 30000,00 | 30000,00 | 30000,00 | 30000,00 | 30000,00 30 000,00 

теап 2,00 25,00 15,78 0,75 4,41 10,50 31,50 

51а 0,82 10,77 9,10 0,94 4,40 2,87 18,19 

тт 1,00 10,00 0,91 0,00 0,07 6,00 3,00 

тах 3,00 45,00 37,50 10,04 38,02 15,00 60,00 


Наиболее лучшие параметры для модели К-ближайших соседей были определены в результате работы 
5-ти блочной кросс-валидации (таблица 5). 


Таблица 5 
Наилучшие параметры модели А-ММ 
Параметр Количество соседей (К) Размер листа (и) Интервал (р) Функция веса (7) 
по 3 15 2 'Иапсе' 
т” 5 15 2 'Чзапсе" 


* 
Наилучшие параметры модели ЗУК для параметров вязкости Мо (в начале релаксационного процесса) и 
модуля скорости 71* были получены эмпирическим путём (таблица 6). 


Таблица 6 
Наилучшие параметры модели ЗУК 
Параметр Тип ядра Порядок ядра Квадратичный регуляризатор Е 
о ‘тЕ 5 0,3 
т“ ОР 6 0,3 
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Соотношение между реальными и предсказанными значениями для модели ^-ММ по параметрам «Вязкость» 


и «Модуль скорости» показаны 


Прогнозируемые значения 


Прогнозируемые значения 


на рис. 3, 4. 
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Рис. 3. Графики ошибок прогноза А-ММ, «Вязкость» 


14 -- =——=ы=-› 
о ый 
12 _  мивеириретиииытиризнит 


10 . 


Актуальные значения 


Рис. 4. Графики ошибок прогноза А-ММ, «Модуль скорости» 


Соотношение между реальными и предсказанными значениями для модели ЗУЕВ по параметрам «Вязкость» и 
«Модуль скорости» показаны на рис. 5, 6. 
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Рис. 5. Графики ошибок прогноза ЗУВ, «Вязкость» 
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Рис. 6. Графики ошибок прогноза ЗУК, «Модуль скорости» 
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Метрики разработанных моделей А-ближайших соседей и опорных векторов представлены в таблицах 7 и 8 


соответственно. 
Таблица 7 
Метрики разработанных моделей А-ММ 
Параметр МАЕ МЕ ВМ$Е МАРЕ, % Е? тат 2 4254 
о 1,8 6,8 2,6 5,9 1,00 0,98 
т” 0,7 0,8 0,9 6,9 0,99 0,98 
Таблица 8 
Метрики разработанных моделей ЗУК. 
Параметр МАЕ МЕ ВМ$Е МАРЕ, % Е? тат 2 4254 
10 1,67 5,75 1,67 8,92 0,98 0,97 
т” 0,72 1,21 1,1 7,3 0,89 0,87 
Помимо синтетических данных, разработанные модели также апробировалась на реальных 


экспериментальных данных, представленных в работе [13]. Использовались экспериментальные кривые 


релаксации поливинилхлорида для различных температур в интервале от 20° до 60 °С. На рис. 7 маркерами 


отмечены экспериментальные значения напряжений при различных температурах в различные моменты 


времени, 


* 
спрогнозированных моделями величин т” и По. 
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Рис. 7. Результаты апробации модели на экспериментальных данных 


а сплошными линиями показаны кривые релаксации напряжений, построенные на основе 
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Обсуждение и заключение. Из рис. 5 видно, что качество прогнозирования по экспериментальным данным 
достаточно высокое, особенно для температур 30°С, 50°С и 60°С. Для других температур качество 
прогнозирования несколько ниже, что связано с качеством самих экспериментальных кривых. Следовало 
продлить время эксперимента и дождаться выхода кривых на горизонтальную асимптоту. 

В данном исследовании наиболее предпочтительным выступает метод опорных векторов (ЗУМ). Это связано 
с тем, что ЗУМ может обрабатывать данные с большим количеством признаков, что является важным для анализа 
реологических параметров материалов. Кроме того, ЗУМ работает с нелинейными зависимостями между 
признаками, применим для решения задачи регрессии, которая является необходимой для определения 
реологических параметров материалов. 

Однако метод СаВооз{ также может быть эффективным в этой задаче, особенно если в данных есть 
категориальные признаки. Кроме того, Са Воо${ может обрабатывать пропущенные данные, что может быть 
важно для анализа реологических параметров материалов. 

Метод А-ближайших соседей является менее предпочтительным в данной задаче по причине невысокой 
эффективности в обработке большого числа признаков, а также наличия проблем с высокой размерностью данных. 

В ходе исследования было показано, что применение методов машинного обучения позволяет эффективно 
анализировать и обрабатывать большие объемы данных, включающие информацию о характеристиках 
полимеров и их реологических свойствах. Модель, разработанная на основе такого анализа, демонстрирует 
высокую точность прогнозирования реологических параметров поливинилхлорида, что подтверждается 
результатами кросс-валидации и сравнением с экспериментальными данными. 

Одним из ключевых преимуществ данного подхода является возможность автоматизации процесса 
прогнозирования реологических параметров полимеров, что позволяет сократить время и затраты на 
исследование и разработку новых материалов. Кроме этого, модель может быть легко адаптирована для анализа 
других типов полимеров и предсказания их свойств. 

В результате выполнения данной работы была разработана прогнозная модель для оценки реологических 
параметров поливинилхлорида методами искусственного интеллекта на основе данных о его характеристиках и 
реологических свойствах. Модель демонстрирует высокую точность прогнозирования и может быть 
использована для оптимизации процессов производства и разработки новых материалов на основе полимеров. 
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